Localized Collective Excitations in Doped Graphene in Strong Magnetic Fields 
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We consider collective excitations in graphene with filled Landau levels (LL's) in the presence of 
an external potential due to a single charged donor or acceptor A~ impurity. We show that 
localized collective modes split ofi' the magnetoplasmon continuum and, in addition, quasibound 
states are formed within the continuum. A study of the evolution of the strengths and energies of 
magneto-optical transitions is performed for integer filling factors v = 1,2,3, i of the lowest LL. We 
predict impurity absorption peaks above as well as below the cyclotron resonance. We find that the 
single particle electron-hole symmetry of graphene leads to a duality between the spectra of collective 
modes for the and A~ . The duality shows up as a set of the D"*" and A~ magneto-absorption 
peaks having same energies, but active in different circular polarizations. 

PACS numbers: 73.20.Mf, 71.35.Ji, 03.65.Ge 
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I. INTRODUCTION 

Graphene, a novel two-dimensional form of carbon,^ 
displays exciting new physics, distinct from that of the 
two-dimensional electron gas (2DEG). A striking exam- 
ple in the presence of a magnetic field is the anomalous 
integer quantum Hall effect, which has been observed at 
room temperature.^'"^ The treatment of electronic inter- 
actions is also challenging in graphene, due to the struc- 
ture and symmetry of the dispersion relations at the 
two inequivalent Dirac points. When scattering between 
these points is negligible, the chirality of the electrons 
results in a suppression of backscattering.A This together 
with graphene's reduced level of defects and impurities, 
makes it highly efficient at charge transport and a promis- 
ing candidate for use in nanotechnology.-^ However, in 
order for this potential to be realized, we need to bet- 
ter understand the precise nature of defects in graphene. 
Optical measurements are a particularly useful tool for 
probing thisi^'^i'SiiiS 

In this paper, we study the magneto- optical response of 
graphene in a strong perpendicular magnetic field in the 
presence of a low density of charged impurities. Infrared 
studies of Landau level (LL) transitions have reported 
significant departures from the bare (non-interacting) 
cyclotron resonance.-^ Whether this can be attributed 
solely to interaction effects and a possible role of dis- 
order remains unclear. Here we develop a general for- 
malism for studying localized collective modes of mag- 
netoplasmon and spin-wave types and determine their 
optical signatures. We treat the electron-electron (e- 
e) interactions beyond the conventional mean-field/RPA 
level and the electron-impurity interaction exactly in 
the high magnetic field regime. Our results are robust 
for a range of impurity screening strengths. We es- 
tablish the existence of an exact symmetry for collec- 
tive excitations, which should be observable by magneto- 
optical spectroscopyii^i^iiS This symmetry, briefly speak- 



ing, connects the magneto-optical electron-like excita- 
tions of, e.g. a positively charged donor with hole- 
like excitations for a negatively charged acceptor A~ . 
This duality is a consequence of the electron-hole sym- 
metry between single-particle states in the lower and up- 
per cones of grapheneJ^ Furthermore, we establish ex- 
act optical selection rules, which demonstrate that the 
"dual" collective excitations with total angular momenta 
Mz = ±1 are active in two different circular polariza- 
tions and, besides having the same energies, exhibit 
the same oscillator strengths. Therefore, a qualitative 
distinction of graphene from the conventional 2DEGfi^ 
is that there are strong dipole-allowed transitions in both 
circular polarizations sensitive to the charge of impurity. 
We show an example of this symmetry for the lowest LL 
n = in Fig. [TJ Each LL in graphene consists of four 
sublevels, due to spin and valley (pseudospin) splitting. 
We denote by {u) a many-electron ground state corre- 
sponding to the sublevel filling factor of a particular 
LL. For sublevels = 1,2,3 of LL with number n, the 
eigenstates and eigenenergies of excitations localized on 
the with M^, coincide precisely with those with —M^ 
formed at filling factor u — A oi the LL with number — n, 
localized on the . This effect represents what remains 
of the electron-hole symmetry after it has been broken 
by a charged impurity. 



II. THE THEORETICAL APPROACH 

We consider collective excitations in a system of elec- 
trons interacting via a screened Coulomb potential Uee = 
e|rf-r2| [13. 14,1^1 in the presence of an additional ex- 
ternal field V . All the results presented here are for a 
Coulomb impurity, V ~ ±e^/eimp|r|r^ but the method 
is applicable for any non-singular axially symmetric po- 
tential T^(|r|). The parameters e and eimp are the effec- 
tive dielectric constants for e-e interaction and electron- 
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FIG. 1: Magnetoplasmons bound on (a) a charged donor at u = 1 and (b) a charged acceptor A~ at u — S. Energies are 
given relative to hujc in units of Eo and e/eimp = 1 (see text). The spectra exhibit the symmetry Z)+ ^ , Mz <-> — Mz and 
<-> 4 — z^. The hatched area of width 0.75-Eo represents the continuum of extended magnetoplasmons. Quasibound states 
within the continuum are not shown. Insets show four branches of resonantly mixed inter-LL transitions conserving spin and 
pseudospin. 



impurity screeningfH respectively. A composite index 
J\f = {nsa} is used to designate the LL number n, the 
spin s =t,i, and pseudospin a =11", -IJ- projections. Low- 
energy collective excitations correspond to the promotion 
of one electron from one of the uppermost filled levels JV2 
to a higher lying empty level A/i (see Fig. [ij. 

Electrons in graphene follow a linear dispersion rela- 
tion close to the zeroes of energy (Dirac points), which 
occur at two inequivalent points in the Brillouin zone, the 
K and K' valleys. We describe a perpendicular mag- 
netic field B by the symmetric gauge A = x r, consis- 
tent with the axial symmetry of T^(|r|). A single electron 
wavefunction in, e.g. the K valley (pseudospin ff), is a 
four-component spinor 

= a„(s„(?!)|„|_i„i(r),(/)|„|„i(r),0,0)xs . (1) 

Here, n is an integer LL number, 0„m(r) is a 2DEG wave- 
function with oscillator quantum number m ~ 0, 1, . . ., 
Gn = 22(^"-o~^\ s„ = sign(n) (with sq — 0) and Xs 
is the spin part corresponding to two spin projections 
s =1, J,.^^ The wavefunction in the K' valley (pseudospin 
4) is obtained by changing the order of the spinor com- 
ponents. The spinors are the eigenstates of the general- 
ized [orbital plus isospin (sublattice)] angular momentum 
projection j^e = Iz + \<yz [l9t with half-integer eigenval- 
ues jz = |n| — m — i. The single-electron energies are 
given by ej^ = sig-a.{n)hjJc\/\n\ + hiOgSz + huycrz, where 
HlOc = vpyJIehB I c is the cyclotron energy in graphene, 
hws is the Zeeman splitting and hio^ is a possible val- 
ley splittingi^S Using the hole representation for all filled 
levels, CMra d}j^^ and cj^„ d^/m for e^/ < ep, we 



introduce operators of collective excitations as 

00 

Qif^Af^M,^ ^A/'iA^2M,(TOi,TO2)cjv'^,„^djv,„, (2) 

mi ,7712—0 

with expansion coefficients satisfying the condition 

^MA/2M,(TOl,m2) ~ dM,,\ni\-Tni-\n2\+m2- An CXact 

quantum number Mz is an eigenvalue of the total Jz = 
ize -\- jzh, for neutral collective excitations Mz is integer 
and of purely orbital nature. It has a direct geometri- 
cal meaning^^ determining the average positions of the 
electron and the hole relative to the impurity, i.e. 

{J^iAf2Mz\rl^rl\AfiM2Mz) 

= (2[M, + 2(|n2|- |ni|)-|-l]+,5„,.o-5„,,o)4 , (3) 

where the states are defined as Q^j\f m \^) = 
W1N2MZ). 

Considering the total Hamiltonian H = Hq + Uee + V 

matrix elements {JViJV^Mz\H\MiJV2Mz) = H^l^^{Mz), 
we find that the effective Hamiltonian is given by 

00 

771 = 
00 

777 = 

'^'^7Vi777i AA27772 '^J\rim[ "'U^m'^ "W27n2 CNimi ■ 

nil -7772 
777^ ,7772 

Here ej\f — ej^ + Ese{J^) denotes the single-particle LL 
energy renormalized by e-e exchange self-energy correc- 
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tions EsE (A/") [11] . Since Kohn's theorem is not applica- 
ble in graphene, these corrections lead to the renormal- 
ization of the bare cyclotron energy, hujc = fi^c + Shwc 
013. For the n = —^ n = 1 transition (denoted here- 
after by Tgi), ShuJc is due only to exchange interactions 
with the lower cone and Shujc — 0.92 Eq- Here Eq = 
(7r/2)-'^/^e^/e/s is the characteristic energy of Coulomb 
interactions in strong B, Ib = {hc/eBY^'^ . 

Due to the spinor form of the single-particle wave- 
functions, the impurity matrix elements in graphene 
are connected with those in the conventional 
Vnm = {<p7im\V{r)\(l)nm), according to 

Va/-™ = (*A/-m|^(?-)|*A/-m) = (s?iV|„|_i,„ + V|„|„,) . 

(5) 

The two-body interaction in ^ consists of the direct 
electron-hole (e-h) attraction and exchange e-h repulsion, 
i.e. 



^.yN'lm\J\r2m2 , . ,7V2m2 AA^m'j^ /„n 



In electron representation. 



2^2/ 



note that we neglect the intervalley scattering in 
graphene by long-range (Coulomb) potentials. Therefore 
the two-particle graphene matrix elements are given by 



nimi 712^2 



^ni ^ni 



I'l |n2l ™2 
|ni I mi \n2 \ m2 



'"i*"! Iml-lmi |n2| "12 "2*n2^|ni| mi |n2|-lm2 ^ ' 
+ S„iS„25niSn2'^|„i|_lmi |„2|_im2 



where Unlmln2m2 — ('/'njm'j ^rijmj l^ee |0nimi '/'n2ni2) '''-'^^ 

those used in the conventional 2DEG. We compute the 
matrix elements for lowest LL's analytically^^ and those 
for arbitrary LL's numerically using Eq. ([8|). 

In general, an infinite number of excitations ^ with 
the same are mixed by the e-e interactions. However, 
those with different single-particle cyclotron energies are 
only weakly (~ Eq/HlOc) mixed in strong magnetic fields 
in graphene and can be neglected. Let us suppose all 
LL's with n < are completely filled, all LL's with n > 
are empty, and the four LL's with n = become succes- 
sively completely filled. We designate the corresponding 
filling factors as = 1,2,3,4. For each i^, there are six- 
teen possible inter-LL excitations involving the n = LL 
as an initial or final state and which have single particle 
energies of magnitude ~ KjJc- Here we concentrate on 
the four excitations for which no spin or pseudospin-flip 
occurs (see the insets in Fig. [1]), as these are the only 
excitations which are optically dipole active. These have 
the same single particle energy hLjc and, therefore, are 
strongly (resonantly) mixed. 



Let us discuss some general features of our approach. 
In the absence of an external potential all magnetoplas- 
mon states are extended and the corresponding Hamilto- 
nian matrix is infinite, see Eq.(|3]). The magnetoplasmons 
can be labeled by a continuous quasimomentum K and 
their eigenenergies fill a band of width ^ Sq.— >^'^ In 
the presence of an impurity, however, some states become 
localized. Importantly, the basis states ^ are localized 
two-particle orbitals whose distances from the impurity 
increase ^ (2m)^/^ZBfi^ Hence, for localized excitations 
the scheme is convergent so that the basis can be trun- 
cated. We include the first TV = 50 basis states for each 
excitation Q^j^-^j^^m with the total matrix size being AN 
for four strongly mixed excitations. The achieved accu- 
racy in eigenergies of bound states is better than 0.1%. 



III. RESULTS AND DISCUSSION 

Figure [ija) shows for v = \ four low-energy branches 
of magnetoplasmons bound on the £)+ for > 0; two 
of these branches are degenerate. For large positive Mz , 
the hole is on average much farther from the impurity 
than the electron, see Eq. Therefore, the e^-D"*" at- 
traction dominates over the h^-D'^ and e-h interactions. 
Generally, for an excitation with the electron in the n*'' 
LL, we find branches with asymptotic Mz 3> 1 ener- 
gies equal to Vnm , when counted from the corresponding 
renormalized cyclotron energy. As an example, notice the 
three branches approaching energy — 0.25i?o and the sin- 
gle branch approaching zero energy in Fig. [TJa). These 
originate, respectively, from the three n = — 1— !-n = 
transitions (denoted hereafter as T-io) and from the sin- 
gle Tqi transition for v = 1. Similar asymptotic behavior 
can be seen for other filling factors. The high-energy (i.e., 
above the band) magnetoplasmons develop for Mz < 0, 
when the hole is on average closer to the £)+ than the 
electron. Such unusual excited states are bound in 2D 
because of the confining effect of Bii^ Fig. [T]^b) shows 
the spectra for j/ = 3 for the A~ which is a "mirror re- 
flection" of Fig.[l](a) because of the aforementioned sym- 
metry. Generally, due to the symmetry, results for the 
at J/ = 1, 2, 3 can be obtained from those for the £)+ 
by changing Mz —Mz and ly ^ 4 — i^. For all filling 
factors, the spectra of bound states are qualitatively sim- 
ilar to those shown in Fig. [TJ For v — A, the states can 
be classified according to the total spin and pseudospin, 
so the states are either spin and pseudospin singlets or 
triplets.^- Only spin and pseudospin singlets are opti- 
cally active. 

Let us consider the magneto-optical response in 
graphene.—"^ In the electric dipole approximation, the 
interaction of electrons with light of frequency iv and left 
(+) and right (— ) circular polarizations is described by 
the Hamiltonian 



evpSe 



ILOC 



(9) 
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FIG. 2: Evolution with filling factor v of energies and optical strengths of magnetoplasmons bound on the with (a) — 1 
active in the ct"*" polarization and (b) with = — 1 active in the a~ polarization. The optically active states are indicated by 
circles with sizes ~ the strongest branches are shown by solid circles (•). The diamonds represent optically dark states. 

The dotted lines are guides to the eye. Inset: Dipole strength \d'Z^ vs energy for v = 2. The spectra were convoluted with a 
Gaussian of width 0.03i5o. The arrow indicates an impurity-related feature below KI)c (below energy zero in the Figure). 



where E is the electric field amplitude and (j± = ioy 
are the PauU matrices acting in the space of the graphene 
crystal sublattices. The exact optical selection rules for 
the collective excitations that follow are: only those with 
no spin or pseudospin flips and with = ±1 and 
1^1 1 ^ 1^2! = il are optically active in the two circular 
polarizations cr*. We quantify the rate of microwave ab- 
sorption in the polarization by calculating the dipole 
transition matrix elements — \ (jvlz — ±l|(5i?±|i^)p 

to final states of magnetoplasmons obtained by numerical 
diagonalization. 

Figure [2] shows optical properties of states bound on 
D"*" for the and a~ polarizations. The results for the 
A~ at = 1,2,3 can be obtained from those reported 
here with the change 1/ ^ 4 — 1/ and ^ . Two 
types of localized states can be optically observed: (i) 
truly bound states, which are split off the continuum 
and have normalizable wavefunctions, (ii) quasibound 
states within the continuum, which have high probabil- 
ity amplitudes on the impurity and long-range oscillating 
tails. The latter may exhibit asymmetric Fano-type op- 
tical signatures 

For both polarizations, the upper branch originates 
mostly from the Tqi transitions with some small (zero 
at J/ = 4) admixture of the T_io. With increasing i', 
the number of Tqi transitions increases, which leads to 
the enhanced contribution of the repulsive e-h exchange 
interactions. This explains the blue shift of the upper 
branch to higher energies with increasing ly. Also, its 
optical strength increases (Fig. [2^) while \d~\'^ de- 

creases (Fig.[5]3). This is explained by the fact that in the 
(7+ (c^) polarization only the electron-like Tqi (hole-like 
T_io) transitions are optically active. 

Screening was shown to be relevant in graphene, al- 
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FIG. 3: Dependence of energies and optical strengths of mag- 
netoplasmons on ratio of effective dielectric constants e/eimp 
for the D"*" in polarization for v — 3. Symbols are as in 
Fig. El 



though the situation is not yet fully understood, partic- 
ularly for a strong magnetic field. Coulomb impurities 
with charge Z = 1 belong to the subcritical regime and 
hence screening effects due to the substrate and the elec- 
tron system in graphene can be modelled via an effective 
charge Zcff < Z^'' or by an effective dielectric constant 
£imp- However the value relative to the e-e screening 
constant, e/eimp is unknown for graphene f^^i^^ We use 
^/eimp = 1 in Figs. [T] and [2] and show in Fig. [3] how 
the energies of the excitations for, e.g. the v = 3 case 
of Fig. mja), are modified by e/eimp < 1- As this ratio 
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increases, the branches of optically active bound states 
are pushed away from the band of extended magneto- 
plasmons. 

IV. CONCLUSIONS 

In conclusion, we established the spectra and the sym- 
metries of collective excitations bound on charged impu- 
rities in graphene in magnetic fields. Our single-impurity 
theory is applicable to samples with finite impurity den- 
sity nimp < I/tt/^, i.e. when the mean separation be- 
tween impurities exceeds the size of bound magneto- 
plasmons. The intensity of impurity peaks will then 
be ^ "-impMii^P- Recent progress in fabrication of 
large graphene films^- with sizes exceeding the wave- 
length A ~ l-Kcj^ujc (> 60 — 100 /zm) opens the way 
for detailed studies of the effects predicted here. Po- 
larization resolved magneto-optical spectroscopy and cy- 



clotron resonance detection using the photoconductive 
response may be very effective experimental probes. Our 
results demonstrate the breaking of particle-hole symme- 
try in a sample with predominantly positive or negative 
impurities, which may partly explain the observed asym- 
metry in LL transitionsi^ The developed method can be 
extended for defects with short-range potentials, results 
will be published elsewhere<^ 
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